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CLAIMS 

1. Method for estimating the seismic ill\amination 
fold \{x,p) in the migrated 3D domain at at least one 
image point for at least one dip of vector p 
wherein the illumination fold I(x.p;s.f) for each (source s, 
receiver F) pair in the seismic survey is estimated, by 
applying the following steps: 

- determination of the reflection travel time 
tr ( 3c^ ^p) is ,r) from the source s to the specular 
reflection point on the plane reflector passing 
through the image point x and perpendicular to the dip 
vector p and then returning to the reflector r ; 

starting from the diffraction travel time td(3c;J,r) 
from the source s to the said image point x and then 
returning to the reflector r ; 

- incrementing the said illumination fold I{x,p;s,?) 
related to the said (source s , receiver f) pair as a 
function of the difference between the diffraction travel 
time td{x;^,r) and ■ the reflection travel time 
tr ( ( ^ ) / s , ; ) . 

2. Method according to claim 1, comprising the step 
of summating each of the said illumination folds I0t.p:i.r) 
related to a (source s, receiver F) pair so as to 
determine the total illumination fold I(3c,^) = ^i(x.p;i,?) . 

s,p 

3. Method according to one of the preceding claims, 
25 wherein, during the incrementing step, the illumination 

fold I(x,p,J,F) is incremented using an increment function 
i(td, tr; ?,F) according to I(x,p) = l(x,p) + i(td,tr; s,r) , the 
said increment function taking account of the difference 
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between the diffraction travel time tci(3c; s, r ) and the 
reflection travel time, tr ( jcj { p ) ; s, r). 

.4. Method according to claim 3, wherein the 
increment function i is a function of the seismic wavelet 
5 s (t) . . 

5. Method according to claim 4, wherein the 
increment function i is expressed as a function of the 
derivative of the seismic wavelet s (t) according to: 
i(td, tr; s, r) = s(td(x; s, r) - trix^i p ) ; s ,r ) 
10 6. Method according to claim 4, wherein the 

increment function i is expressed as a function of the 
derivative s (t) of the seismic wavelet s (t) with respect 
to time according to: 

i(td/ tr; s, r) = td(x; s, r) - t^ix^ip ) ; s ,r ) 
15 7. Method according to any one of claims 3 to 6, in 

which an a priori correction w(jc,J,F) of the illumination 
fold is taken into account by migration, comprising the 
step of incrementing the illiamination fold l{x , p ;s ,r) 
related to a (source s, receiver F) pair by 
20 i(td,te;^ ,r ) .w(3c;s , r ) . 

8 . Method according to any one of the preceding 
claims, wherein the determination step includes the 
second order Taylor series development of the diffraction 
travel time td(x;J,r) around the image point x: 

25 td(^;J,r )=td{x; J,r ) + (V«td(3c;?,F))'. (3c,-jc)+- (5c,- 
^)'-Ax,xtd(X; J, r ) . (Xr-x) 

9. Method according to claim 8, wherein the specular 

reflection point (P ) is determined along the length of 
the said reflector such that the. diffraction travel time 

30 at the said specular reflection point (P ) is 
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Stationary^ according to the equation: 

p'^ A( Vxtd(x;?,r )+Ax,«t<i(X;^,r ) . ( 3c^ ( ^ ) - 3c ) ) = 0 . 

10. Method according ■ to . any one of claims 8 or 9, 
wherein the specular reflection point and the 

5 reflection travel " time tr(Xr(p);s,r) (are determined 
according to the following expressions: 

^r(^) = X'M.F'^. b 

trixrip) fS ,r ) = td(3c ;?,r .fe"^. F'^. b 

2 

10 

where : 

- M is a (3x2) matrix described by two vectors 
extending along the length of the reflector, and 
therefore perpendicular to the dip vector p ; 

15 -feisa(2xl) vector of first order derivatives 

of the diffraction travel time along the reflection 
plane: b = m'^. ( V^td) ; 
t • ~ F is a (2 X 2) matrix of second order derivatives 

of the diffraction travel time along the reflection 

20 plane: F = M^. (A .^^^td) .M. 

11 . Method according to any one of the preceding 
claims, wherein the determination step uses isochronic 
migration maps ta{x;s,r) specified for each (source 
receiver F ) pair involved in the migration at each image 

25 point X in the migrated 3D domain. 

12. Method according to any one of the preceding 
claims, wherein the seismic illumination fold l{x,p) in 
the migrated 3D domain is estimated during the Kirchhoff 
summation migration of seismic data recorded during the 

30 3D seismic prospecting. 
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13. Method for correction of seismic data amplitudes 
recorded during ' 3D seismic prospecting in order to 
compensate for the effect of non-uniform illumination of 
sub-soil reflectors, comprising the steps of: 
5 - estimating the illumination fold l(x,p) using the 

method according to any one of claims 1 to 12, 

- using the inverse I'^{x,p) of the said ratio as a 
weighting factor to be applied to each of the said 
seismic data amplitudes. 

10 ■ 14. Method for selection of an acquisition geometry 

among a plurality of acquisition geometries as a function 
of the target of 3D seismic prospecting, comprising the 
steps of: 

- determining the illumination fold l{x,p) by the 
15 method according to any one of claims 1 to 12, for each 

of the acquisition geometries considered, 

- selecting the acquisition geometry providing the 
optimum illumination fold as a function of the target. 



